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We present first measurements on a superconducting three-dimensional, partly chaotic microwave 
billiard shaped like a small deformed cup. We analyze the statistical properties of the measured 
spectrum in terms of several methods originally derived for quantum systems like eigenvalue statis- 
tics and periodic orbits and obtain according to a model of Berry and Robnik a mixing parameter of 
about 25%. In numerical simulations of the classical motion in the cup the degree of chaoticity has 
been estimated. This leads to an invariant chaotic Liouville measure of about 45%. The difference 
between this figure and the mixing parameter is due to the limited accuracy of the statistical analy- 
sis, caused by both, the fairly small number of 286 resonances and the rather poor desymmetrization 
of the microwave cavity. Concerning the periodic orbits of the classical system we present a com- 
parison with the length spectrum of the resonator and introduce a new bouncing ball formula for 
electromagnetic billiards. 

PACS number(s): 05.45. +b, 41.20.Bt, 84.40.Cb 



I. INTRODUCTION 

In the last few decades the theoretical investigation of 
two-dimensional Euclidian and Ricmannian geometries, 
so-called billiards, has led to a very fruitful new disci- 
pline in non-linear physics |l]-|j|. Due to the conserved 
energy of the ideal particle propagating inside the bil- 
liard's boundaries with specular reflections on the walls, 
the plain billiard belongs to the class of Hamiltoman sys- 
tems with the lowest degree of freedom in which chaos can 
occur and this does only depend on the given boundary 
shape. Because of their simplicity two-dimensional bil- 
liards are in particular adequate to study the behaviour 
of the particle in the corresponding quantum regime 
where spectral properties are completely described by the 
stationary Schrodinger-equation 
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inside the domain Q with Dirichlet boundary conditions 
on the walls 
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In this context the investigation of " Quantum Chaos" has 
become one of the most fascinating goals of theoretical 
physics at the end of this century . 

About five years ago experimentalists have even found 
very effective techniques to simulate the quantum billiard 



problem with the help of macroscopic devices. Due to the 
equivalence of the stationary Schrodinger-equation and 
the classical Helmholtz-equation in two dimensions one 
is able to model the billiard by a similarly shaped electro- 
magnetic cavity p|-pT[. In former publications we have 
demonstrated the high accuracy of large ensembles of 
measured eigenvalues as well as of resonance shapes and 
the attached widths in superconducting cavities formed 
like desymmetrized Bunimovich and truncated Hyper- 
bola billiards ||-|l| . 

Especially the statistical analysis of measured (or nu- 
merically simulated) eigenvalue sequences confirms that a 
distinction between classical chaotic and regular systems 
from the quantum point of view is only possible in the 
semiclassical regime (formally spoken for Ti — > 0) where 
the particle's de Broglie wavelength is sensitive to details 
of the billiard's borderline. One of the most surprising 
and primarily empirical results of these investigations is 
the fact that quantum spectra of classical chaotic Hamil- 
tonian systems can be described in a universal manner 
which only depends on the global symmetry of the under- 
lying dynamics: The quantum pendants of time-reversal 
invariant classically chaotic systems typically reveal spec- 
tral structures which are reproduced excellently by sta- 
tistical properties of the Gaussian Orthogonal Ensemble 
(GOE) of Random Matrix Theory On the other 

hand classically regular systems usually lead to spectral 
fluctuations on the wave dynamical side according to un- 
correlated Poissonlike distributed random numbers fllyl. 
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In this article we present first investigations which were 
performed on a three-dimensional superconducting bil- 
liard. Due to the polarization properties of the electro- 
magnetic fields E and B inside the cavity the full vecto- 
rial Helmholtz-equations pfj 



(A + e^-j)E(f) = 



(A + eu- J )B(r) = 



(3) 
(4) 



have to be used with corresponding boundary conditions 



E\\(r)\ d g=0 and B x (f)\ dg = 



(5) 



on the walls which are assumed to be ideally conduct- 
ing. Of course the analogy with the corresponding scalar 
Schrodinger-equation in the same geometry is fully lost. 
Instead of talking of the semiclassical limit we have to de- 
scribe the classical electromagnetic billiard in this region 
in terms of ray-optical characteristics, where features of 
the periodic orbits inside the geometry dominate the cor- 
responding wave optical side. 

Three-dimensional systems have so far only scarcely 
been investigated experimentally. To our knowledge first 
experiments with electromagnetic waves in cavities simu- 
lating acoustic wave phenomena in rooms were performed 
by Schroder pi]] . Acoustic model statistics in metal 
blocks have also been investigated [ p2|j2^1 . Very recently 
statistical properties of eigenfrequency distributions in 
asymmetrically shaped microwave cavities have been re- 
ported 24 1. Theoretically quantum effects as well as 
their electromagnetic counterparts in three-dimensional 
systems were treated in Refs. |27| . 



II. EXPERIMENT 

We have investigated a small deformed 3D-cavity made 
from an open deep drawn Niobium cup (a gift, which we 
had received from the director of the Continuous Electron 
Beam Accelerator Facility, CEBAF, Hermann Grunder) 
with a welded lid to cover it. Before welding the two 
parts together, the cup was deformed at its open end 
to destroy its rotational symmetry. Its measures are 
given in Fig.l. The total volume of the resonator was 
determined by filling it with water and measuring the 
additional weight, yielding a volume of (122.5 ± 1.9)cm 3 . 
The shape of the cup can be approximated by cutting a 
three-axial ellipsoid twice perpendicular to its main axis 
(see also Fig.l) and then deforming it slightly. The mi- 
crowaves were transmitted into and out of the cavity by 
two small antennas in the lid. The whole resonator was 
cooled down in Helium-atmosphere at a temperature of 
2K and a pressure of 38mbar in one of the cryostats of the 
superconducting Darmstadt electron linear accelerator S- 
DALINAC |Q together with the accelerating structures. 
With this setup we have measured microwave spectra in 



transmission (antenna 1 for input and antenna 2 for out- 
put) as well as in reflection (same antenna for excitation 
and detection) using a Hewlett Packard network analyzer 
(model HP8510B) in a frequency range between and 
20GHz, respectively. 

Spectra were taken in 10kHz steps and Fig. 2 shows an 
extraction of the transmission spectrum between 15 and 
20GHz where typical Q- values of up to 10 5 and signal-to- 
noise ratios of up to 50dB were obtained. By comparing 
the three measured spectra 286 resonances could be con- 
sistently identified which form the base of all following 
investigations. 

For an analysis of this set however one has to take into 
account the following point: Since the geometry of the 
cup is very close to one which possesses two symmetry 
planes (xz— and yz— plane, see Fig.l) the given set of 
resonances in the analysis is always compared to a super- 
position of four independent subspectra. 



III. RESULTS AND DISCUSSION 

A. Density of eigenmodes 

In order to derive meaningful statistical measures for 
the given eigenvalue sequence it is first necessary to ex- 
tract the smooth part of the resonator's density of eigen- 
modes which is given by the generalized electromagnetic 
Weyl-formula (29|-|32| 
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(f) = ^\G\f + const. 



(0) 



where \G\ denotes the volume of the cavity and / the up- 
per frequency limit of the given spectrum. The constant 
term contains contributions of the surface's curvature as 
well as of the edges of the cavity. The total density of 
eigenmodes contains in addition a fluctuating part 



p(f) = P s 



olh 



(f)+pf^f) = J2 S(f-fi), (7) 



where ft denotes the eigenfrequencies of the resonances. 

It is very instructive to compare Eq.(^|) with 
the corresponding expression for the three-dimensional 
Schrodinger-problem in the same geometry, i.e. 
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\dQ\f + const. 
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Due to the already mentioned polarization features of the 
modes, the leading term in the electromagnetic formula, 
Eq. (^|) , is twice the corresponding term in the scalar prob- 
lem, Eq.(|J), because of two transverse directions of po- 
larization relative to the axis of propagation. In addition 
these two polarizations which are known as TM- and TE- 
modes in certain geometries provide linear terms in the 
smoothed eigenmodc density of the same magnitude but 
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with different sign. Thus due to cancellation there is no 
linear contribution in Eq.(||), whereas this term survives 
in the scalar case and is proportional to the cavity's sur- 
face \dQ\. Even if a clear separation in the attached TM- 
and TE-modes is not possible in cases of arbitrary geome- 
tries, like in the case of our cup, this linear term vanishes 
for all piecewise smooth boundaries |3^j3^] . The constant 
terms in both Eqs. @j and (||) are of the same origin, 
i.e. the curvature of the surface and the edge angles of 
the cavity. 

To determine the spectral fluctuations, the smooth 
part of the eigenmode density in the measured spec- 
trum had to be eliminated. For this we constructed from 
Eq.(^) the staircase function 

N(f) = f dfp(f') (9) 
Jo 

=E (/-/^)=E 1 +E5 

f>fi f=h 

and obtained its fluctuating part 

N fluc {f) = N(f) - N smooth (f) . (10) 

Since in the case of our billiard there is no analytical 
form for the edge contribution of the constant term of 
Eq.(^|) we have fitted a third order polynomial (without 
the quadratic term) to the experimental staircase func- 
tion, i.e. 

N smooth (f) = Vf f + V 2 - f + V 3 . (11) 

In Fig. 3 the remaining fluctuating part of the staircase 
function can be seen to oscillate around zero as ex- 
pected, the fitted constant V\ corresponds to a volume of 
(119.3 ± 0.7) cm 3 which is very close to the correct value 
of (122.5 ± 1.9)cm 3 . This difference in the leading term 
corresponds to an uncertainty of 4 resonances in the mea- 
sured total spectrum (neglecting the linear and constant 
term in Eq.(|ll])). 



B. Nearest neighbour spacing distribution 

In order to perform a statistical analysis of the 
given eigenvalue sequence independently from the special 
size of the resonator, the spectrum was first unfolded, 
i.e. from the measured sequence of eigenfrequencies 
{fi, /2, •••,/>, fi+i, ■ ■ ■} the spacing s 4 = (f i+1 - /;) /a 
between adjacent eigenmodes was obtained by calculat- 
ing the local average s from Eq.(|ll]). The proper nor- 
malization of the measured spacings of eigenmodes then 
yielded the desired nearest neighbour distribution P(s), 
i.e. the probability for a certain spacing s. 

By comparing P(s) to theoretical expressions one has 
to take into account, however, that the mechanically only 



weakly deformed cup (Fig.l) can to a good approxima- 
tion be constructed of four similar quartercups which sep- 
arately possess the full geometrical information of the ob- 
ject, in which case the measured spectrum would be a su- 
perposition of four independent symmetry classes which 
are obtained by permuting the boundary conditions on 
the cutting planes of the quartercup from electric to mag- 
netic, respectively. Hence one has to be aware of this 
point in the following statistical analysis of eigenmode 
spacings. 

Furthermore, to obtain a quantitative criterion con- 
cerning the degree of chaoticity in the system the spec- 
trum was analyzed in terms of statistical measures from 
a model of Berry and Robnik [^3) which interpolates 
between the two limiting cases of pure Poissonian and 
pure GOE behaviour for classical regular and chaotic sys- 
tems, respectively. The final essence of this model is a 
mixing-parameter q which is directly related to the rel- 
ative chaotic part of the invariant Liouville measure of 
the underlying classical phase space in which the motion 
takes place. According to the herein embedded one-to- 
one connection between classical phase space and eigen- 
mode density for the two different regions of regular and 
chaotic motion, a comparison becomes meaningful only 
for the highly excited domain of the spectrum which is 
not sufficiently covered by the frequency range up to 
20GHz investigated here. 

In the upper part of Fig. 4 the result for the near- 
est neighbour spacing distribution P(s) which describes, 
as already mentioned, short range correlations between 
neighbouring unfolded levels is shown in form of a his- 
togram. With respect to the superposition of four inde- 
pendent symmetry classes, as noted above, the model of 
Berry and Robnik makes the following ansatz for an in- 
terpolation between pure Poissonian and pure Gaussian 
characteristics 

P M M = e- (1 - 9)s (l6(^) 2 er/c 4 (^| S ) + (12) 
L 4 4 2 V 4 7 J 

CXp(-J(f) 2 S V/c3(^f S ) + 

12(|)^exp(-2j(|)V)e,/c 2 (^| S )) . 

For inspection the limiting curves of a superposition of 
four independent pure Poissonians (q = 0), which again 
yields one single Poissonian, as well as of four indepen- 
dent GOEs (q = 1) are also represented in the figure. 
Note that the superposition results in a clear loss of lin- 
ear level repulsion for small spacings in the case of pure 
GOE statistics. Furthermore there is not much distinc- 
tion between Poissonian and GOE statistics, anyhow. 

It is obvious that P(s) does not allow to determine the 
mixing parameter q with reasonable significance, since 
the fluctuations in the data are larger than the difference 
between the given smooth curves. To be free of effects 
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due to the binning of P(s) we have also calculated the 
cumulative nearest neighbour spacing distribution 



1(a) = / P(s')ds' 



(13) 



the result is presented in Fig. 4, lower part, together 
with the curves for the pure distributions, Poisson and 
4xGOE. In this case the measured spacings are very close 
to Poissonian behaviour over a large range of s. A fit of 
J(s) to the data yields q — 0.16±q' 2 |, where the uncer- 
tainty was determined from analyzing different subsets 
of the measured data. Because of this rather poor sensi- 
tivity other statistical measures like the number variance 
and spectral rigidity had to be applied to the data to 
determine q. 



C. Number variance and spectral rigidity 

In order to check for long range correlations between 
the measured levels we have calculated S 2 — as well 
as A3— statistics, two measures originally introduced by 
Dyson and Mehta [ ^6|Jl^1 for studies in equivalent fluctu- 
ations of nuclear spectra. In this case one is interested 
in spectral correlations on a scale which contains L mean 
level spacings. Here 

£ 2 (L) := <(n(L)- < n(L) >) 2 ) =< n 2 (L) > -L 2 (14) 

represents the averaged variance of a number n(L) of lev- 
els belonging to an interval of length L on the unfolded 
axis with mean < n(L) > = L. The quantity A^(L) is 
a smoothed and rescaled version of £ 2 (£) and can be 
calculated from 



A 3 (L) = ^j Q dr(L 3 -2L 2 r + r 3 )E 2 (r) . (15) 

Figure 5 shows the experimental results for these two 
measures as well as the fitted curve according to Berry 
and Robnik [1331] deduced from 



T,jJL,q) 



■^Y? GOE (qj) 



(16) 



where T? Poisso n 



and ^qoe are the distributions for the 
pure cases and one symmetry class. An equivalent for- 
mula is valid for Aa(L). As in the case of the nearest 
neighbour spacing distribution the superposition of four 
pure Poissonians yields one single Poissonian which can 
be directly seen from the analytical expression 



AL) = L 



(17) 



Inspecting Fig. 5 two results have to be noted: Firstly, the 
experimental number variance E 2 lies clearly between the 
limiting curves of pure Poissonian and Gaussian charac- 
teristics. A fit of the expression (pi) to the data yields a 



within the given error range with to the g-value derived 
from the cumulative nearest neighbour spacing distribu- 
tion. On the other hand the A3 (L)— curve only corre- 
sponds within the given error range of the data with this 
result, i.e. a deviation from pure Poissonian behaviour 
is basically not visible. Secondly, the number variance 
S 2 (L) clearly displays saturation above a certain value 
L m ax- In fact, Berry showed |34|] that a global and uni- 
versal semiclassical (or ray-optical) behaviour can only be 
expected between L m m = 1 and the value L rnax which is 
determined on one hand by the finite ensemble of reso- 
nances, but above all by the lengths of the shortest pe- 
riodic orbits of the classical system. From Fig. 5 we find 
Lmax ~ 10 for E 2 . Note, that L max of A3 is well ap- 
proximated by four times L max of S 2 (see Eq. (|15|) and 
Pq| ) . Using a modification of Berry's expression jpjfl for 
the electromagnetic case, L max of A 3 can be related to 
an average length lmin of the shortest periodic orbits via 
the expression 
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lmin f ma > 



(18) 



mixing-parameter of £ (£), q — 0.30±q; 2 q, which agrees 



where An is equal to the first term V\ ■ / 3 in Eq. (|Tl|) and 
fmax denotes the upper frequency, i.e. f m ax — 20GHz. 
The result is l m i n ~ 0.17m in fair agreement with a value 
deduced independently in Sect.E below. 

In addition to £ 2 (£) and As(L) we have also calcu- 
lated another long range statistics: the two-level form 
factor 62 (i) in the corresponding time domain |35 36 
which leads to a so-called "autocorrelation hole" in the 
case of a GOE-like sequence. It turns out that the fluc- 
tuations in 62 due to the small number of resonances are 
to large for a proper conclusion and even in the pure 
GOE-case the hole is not very pronounced due to the 
superposition of symmetry classes. 

To summarize the statistical investigation, only the 
cumulative nearest neighbour spacing distribution J(s) 
and the number variance S 2 lead to a significant devi- 
ation from pure Poissonian behaviour for the given set 
of resonances. Furthermore the extracted mixing pa- 
rameter q represents an upper limit for the chaoticity 
since the assumption of a superposition of four symme- 
try classes in the statistical analysis increases the weight 
of the pure GOE-contribution with respect to the desym- 
metrized case. 



D. Classical surface of section 

In order to obtain an independent estimate for the de- 
gree of chaoticity in the system we have also performed 
a numerical simulation of classical motion inside the 3D- 
cup. Following an idea of Zaslavsky and Strauss |37|| 
we have approximated the resonator's geometry by one 
half of a so-called barrel billiard which can be obtained 
by cutting a three-axial ellipsoid in two different heights 
perpendicular to the longest of its axis, see also Fig.l. 
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Consequently the surface of the barrel can be described 
by the curve 

2 2 2 

+ to + = 1 with z e [-ft, 0] (19) 
a z o z c 2 

and the resonator was modelled with axes a = 
58.5mm/2,& = 67.6mm/2 and cuts at z — and z = 
—ft = —56.1mm. With respect to these parameters the 
relative geometrical difference in the volume between the 
cup resonator and the barrel billiard is about 11%. 

To characterize the classical motion of a particle inside 
this "half of a barrel geometry" a two-dimensional area 



preserving mapping also introduced in |37| was used for 
creating the underlying Poincare surface of section. The 
base of this mapping is simply a pair of a coordinate and 
its conjugated momentum (z,p z ) of the particle during 
each reflection on the wall. Note that the momentum 
|p| is normalized to unity. Figure 6 shows the resulting 
patterns in phase space after 16000 collisions with the 
boundary. Because the full ellipsoid is totally regular ]39 ] 
there exists a class of trajectories (also shown in Fig. 6) 
manifested by regular stripes on the surface of section 
which indicate remaining regularities inside the cup. This 
is due to the fact that these trajectories do not impact on 
the plane at z = — h which exactly is the reason for over- 
focussing features of the geometry and consequently the 
origin of chaos. The relative area of the chaotic part of 
phase space, the so-called "chaotic sea" , was estimated by 
an invariant chaotic Liouville measure of q « 45% using a 
drastically increased number of about 2-10 5 wallcollisions 
to get a higher accuracy. This is bigger than the mixing- 
parameter deduced from I(s) and E 2 , but it has to be 
noted that the applicability of the Berry-Robnik-model 
was originally shown to depend strongly upon reaching 
the semiclassical (or ray-optical) limit which in the two- 
dimensional case is established only far beyond thousands 
of eigenvalues |59| . The present case of only close to three 
hundred eigenvalues at least points to the correct ten- 
dency for the correspondence between the classical Liou- 
ville measure and the wave dynamical mixing-parameter 
extracted from the data. 



E. Periodic orbit theory 

As a final study of semiclassical features of our wave 
dynamical system, we have analyzed the spectrum of 
classical orbit length I in the 3D-cup which is directly 
related to the measured frequency spectrum by a Fourier 
transform of the above mentioned fluctuating part of the 
eigenmode density (Eq.(Q)), 

~P fluc {l) = / P fluc (f) ■ exp(i—lf) df (20) 

(P(f) ~ p sm °° th (f)) ■ exp(i— If) df . 

If rain C 



Here f m in and f ma x denote the borders of the measured 
frequency range, i.e. and 20GHz. Figure 7 shows in 
the upper part the value \p' luc (l)\ 2 in a range of orbit 
lengths I up to 0.5m, and in the lower part some periodic 
orbits obtained from numerical simulations on the barrel 
billiard are presented. As can be seen from the arrows 
below the abscissa several of those orbits correspond well 
to the locations of peaks in the Fourier spectrum of the 
data. The shortest periodic orbit which bounces between 
the bottom and the lid has a length of I = 0.1122m, which 
is close to the independent estimate of Sect.C above. It 
is the so-called shortest bouncing ball orbit in the cup. 



F. 3D-bouncing ball orbit 

This bouncing ball orbit has been analyzed in a more 
quantitative manner. Following pp[ | we have calculated 
an additional term of the smooth part of the stair- 
case function, Eq.([ll"|), attached to the three-dimensional 
bouncing ball orbit propagating periodically between the 
bottom and the lid of the cup. Again the contributions 
for both polarization classes have to be considered sepa- 
rately. The result for the sum of both parts as pointed 
out in the Appendix is given by 



(21) 



ttS 
2ft 2 



£ (X 2 -n 2 )-^ 3 
0<n<X 



-X' 



where X = kh/w = 2ft//co and ft = 56.1mm, the half 
of the bouncing ball orbit's length, i.e. the heigth of 
the cup. The parameter S denotes the area on which 
this orbit exists, i.e. the size of the cup's bottom. The 
upper part of Fig. 7 also shows the result for the remain- 
ing length spectrum after extracting this contribution. 
In fact, the peak belonging to the correct length of the 
bouncing ball orbit does not vanish completely because it 
is due to a superposition of two adjacent periodic orbits, 
the bouncing ball orbit at I — 0.1122m and a stable pe- 
riodic orbit at I = 0.1170m, as seen from the lower part 
of Fig.7. 

In order to verify the new electromagnetic bouncing 
ball formula, Eq. (pi|) , independently from this result we 
also have tested it using a set of w 20000 eigenmodes of 
a regular box. This system is classically integrable, thus 
it is possible to calculate the electromagnetic eigenfre- 
quencies analytically. The box also allows to study the 
different polarizations, TE- and TM-modes, in more de- 
tail and especially their systematic degeneration which 
directly follows from the analytical expression for the 
eigenfrequencies 



fu,v,' 



Co 

2 



(J 2 + (|) 2 + (^) 2 with (22) 



TM : u,v= 1,2,3,... w = 0,1,2,... 

TE : u,v = 0(either u or v),l,2,... w = 1,2,3,. 
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In order to avoid further accidental number theoretical 
degeneracies the lengths of the three edges A, B and C 
were chosen as A = 0.2m, B = jA and C — j 2 A, where 
I/7 = (s/5 — l)/2, the ratio of the golden mean. Follow- 
ing Eq.(^) we have calculated the exact smooth part of 
the staircase function for the given set of eigenfrequen- 
cies. Here the linear term as well as the constant term 
are given analytically without any free parameter. Ex- 
tracting this smooth part from the total staircase of the 
box, Eq.(^|), we obtained the fluctuating part, which is 
shown in the upper part of Fig. 8. The present choice of 
the edges A, B and C yields a very clear oscillation in 
this fluctuating part, which is dominated by one certain 
classical bouncing ball orbit (also given in Fig. 8). The 
length of this orbit, Ibbo, corresponds to the period of the 
observed oscillation A/ = Co/hbo — c$/2A w 750MHz. 
Using the area Sbbo — B-C for this certain orbit, we have 



calculated the contribution which follows from Eq.(21) 



This curve is also represented in the upper part of Fig. 8 
and reproduces the data very well as it can be seen in the 
lower part of the same figure where the bouncing ball con- 
tribution has been extracted from the staircase function, 
i.e. where the modulation of A/ is not present anymore. 
Note that the impressive amplitude of the given modula- 
tion is due to the ratio 2ir ■ Sbbo/llbo = ' 7 3 ~ 6.65 in 
Eq.(pll) which is the largest contribution in comparison 
to the other possible bouncing ball orbits and also much 
larger than the corresponding amplitude for the 3D-cup, 
nS/{2h 2 ) « 0.45. 

To perform a more obvious test of Eq.( |2l|) we again 
have calculated the Fourier spectrum, Eq. (|20|) , of the 
calculated set of eigenfrequencies. The result is shown 
in Fig. 9. The bouncing ball orbit leads to an impressive 
peak at a length 4&o = 0.4m. Using Eq.(21) in the same 



way, i.e. calculating the Fourier transformed of N^„° we 
are able to extract this peak without any remnant as well 
as contributions of its multiples, which indicates that the 
orbit has been described correctly by Eq.(21). Because 
of this result, which was achieved successfully for several 
boxes of different proportions, we are certain that the 
remaining peak in the length spectrum of the 3D-cup in 
Fig. 7 is due to the existence of the stable periodic orbit 
of 0.1170m length. 

As an additional test of statistical measures, we have 
calculated the nearest neighbour spacing distribution, 
the E 2 — and the A3— statistics for both polarizations of 
the box separately. This yields pure Poissonian charac- 
teristics accross the whole range of spacings as expected 
for the given regular system. As an example, the result 
for the nearest neighbour spacing distribution using the 
first 9978 TE- modes is presented in Fig. 10. 

Return to the 3D-cup: For completion we have cal- 
culated the statistical measures of the short and long 
range correlations after extraction of the bouncing ball 
contribution, i.e. we have reunfolded the experimental 
spectrum in adding Eq. (pl|) to the standard Weylian, 
Eq. (|i"l|) , and repeated our statistical analysis. The ef- 
fect is rather weak, as can be expected because of the 



small amplitude of Eq.(|2i|) in the case of the cup. The 
procedure leads to a very slight correction of the E 2 — 
and the A 3 — curves towards the pure GOE-limit but its 
influence on the mixing parameter is much smaller than 
the quoted uncertainties. 



IV. CONCLUSIONS 



In summary we have applied several original quan- 
tum methods like eigenvalue statistics and periodic or- 
bits in the given case of a purely classical, vectorial 
wave phenomenon. Although the spectrum was checked 
for completeness using the correct electromagnetic Weyl- 
formula, Eq.(ll), the fluctuations around this mean be- 
haviour are not very sensitive for the details of the bound- 
ary in the present case of only 286 resonances. On the 
other hand the system was analyzed in terms of periodic 
orbit theory. As the length spectrum of the resonator 
leads to more detailed conclusions (estimations for the 
saturation of E 2 and A3 via the shortest periodic orbits, 
identification of classicals orbit lengths and extraction 
of the first bouncing ball orbit) this analysis reveals as 
a helpful tool for investigations in the near scmiclassi- 
cal limit. Comparing the classical degree of chaoticity 
(<lci ~ 0.45) with the corresponding mixing parameter 
deduced from spectral statistics according to the model 
of Berry and Robnik, only the cumulative nearest neigh- 
bour spacing distribution I(s) and the number variance 
E 2 significantly deviate from the pure Poissonian be- 
haviour and show the correct tendency (q > 0). The 
uncertainties of the analysis are not due to fundamental 
difficulties but rather a consequence of the poor desym- 
metrization and the small size of the cup, yielding only 
286 resonances in the accessible frequency range. There- 
fore, in order to study the "three-dimensional Helmholtz- 
Chaos" in a more quantitative and proper way we have 
already started measurements on a precisely manufac- 
tured, fully desymmetrized 3D-Sinai-Billiard realized by 
1/48 of a cube with a centered sphere 26 1. Besides test- 
ing the usual statistical measures we will especially inves- 
tigate the ray-optical generalization of Gutzwiller's quan- 
tum trace formula following the very advanced ideas of 
Balian and Duplantier |32|. 
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APPENDIX: THE BOUNCING BALL 
STAIRCASE FUNCTION IN D DIMENSIONS 

First, we construct the general contribution of a bounc- 
ing ball orbit in a Z?-dimensional cavity (D > 1) to the 
spectral staircase function of non-relativistic quantum 
mechanics, N^°(k) with k = 2irf/co, and secondly we 
apply these general results to the electromagnetic case 
in D = 3 dimensions. As usual the spectral staircase 
function is given as the integral over the spectral density, 
Pd°(E), 



N%°(k)= / dkp b £°(k) 
Jo 

ak 



(Al) 



k dk^p%°(E) 
m 



where m is the mass of the non-relativistic particle inside 
the cavity. The spectral density in turn is given by the 
imaginary part of the trace of the non-relativistic Green's 
function inside the cavity: 



1. 



Pd°(E) = - lim -ImTr(G(E + ie)) 



(A2) 



1 s 1 - ^^ 

— lim — Im > t— ; - - x 

e — >o 7r ^ 2n ) L> - 1 



Here I is half of the total length of the bouncing ball 

orbit and sf 1 ' is the size of the smallest of the two 
parallel D — 1 dimensional "surfaces" between which the 
orbit bounces. The integer n labels the modes along the 
bouncing ball orbit with n > or n > for Dirichlet 
or Neumann boundary conditions on the two surfaces, 
respectively. The momentum component along the orbit 
is nir/l. The D — 1 dimensional momentum integration 
perpendicular to the bouncing ball orbit is the free unre- 
stricted one. The factors in front of the integral are the 
phase space normalization factors in D — 1 dimensions. 



After the insertion of Eq.(A2) into Eq.(Al), the spectral 
staircase function reads 



r k* a(D-l) 

N%°(k)= I d(k 2 ' ■ 



tE 



{2-kY 



dk ± k°- 2 8{k 2 , 



dfll ^ x (A3) 



2 n 2 



where /dfij£ -1) = 2tt^ d - 1 ^ 2 /T({D - l)/2) is the 
angular integral in D — 1 dimensions. In writ- 
ing down Eq.(A3) we used the following relation: 



lim e ^o Ith(x — xq — ie)/Tt — 5(x — xq). Now the integra- 
tions in Eq.([A3]) become simple because of the delta func- 
tion. The final result for the bouncing ball orbit of a non- 
relativisitic particle inside a cavity in D > 1 dimensions 
is therefore 



N b D bo (k) = 



S 



(D-l) 



T{^)(2^)^ 

D-3, | 2 E 



(A4) 



2 2 \ CD-3)/2+l 



P 



This formula applied to two dimensions reads 



L 



E 



k 2 - 



(A5) 



where L = S± is the length perpendicular to the bounc- 
ing ball orbit. It agrees of course with the result of Ref. 
p(| . The three-dimensional expression has the form 



N bbo 



(k) 



47T 



E 



k 2 - 



r- 



(A6) 



(2) 

where S = S ± denotes the plane perpendicular to the 
bouncing ball orbit. The bouncing ball contribution in 
three dimensions to the fluctuating part o f th e spectral 
staircase function is now given by formula (A6) minus its 
contribution to the Weyl or smooth part of the spectral 
staircase function which has to be subtracted in order 
to a void double counting. The smooth part follows from 
Eq. ( |A6| ) via the application of the Euler-MacLaurin for 
mu\t^F(0) + F(l) + F(2) + ■■■ F(N - 1) + ^F(N) = 

J N dxF(x) + ... 
dots correspond 
Thus we have 



) + $F(N) 
where F(x) cx (k 2 — x 2 n 2 /l 2 ) and the 
to the fluctuating contribution here. 



Ar&fro, smooth / ; \ 

N 3 (k) 




(+)2 ; 



Slk 3 - Sk? 
6tt 2 (+) 8tt 



(A7) 



where the last term takes into account the factor 1/2 in 
front of the first term, -F(O), in the Euler-MacLaurin for- 
mula and the extra mode in the Neumann case. We do 
not get an extra contribution from the upper boundary, 
as F(N) vanishes in the average in our case. In summary, 
in non-relativistic quantum mechanics, the bouncing ball 
contribution to the fluctuating part of the spectral stair- 
case function in a three-dimensional cavity reads 



N. 



bbo, flue 



,5' 



w = ^ E (* 2 - 



Slk 3 + Sk 2 
(A8) 
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where the upper signs/inequalities apply to Dirichlet 
boundary conditions on the two surfaces between which 
the orbit bounces, whereas the lower signs/inequalities 
refer to the Neumann case. In the electromagnetic case 
these two contributions correspond to the magnetic and 
electric bouncing ball modes, respectively, see Eq.(||), 
which decouple for modes along the bouncing ball or- 
bit. Therefore we can just sum both terms to get the 
final expression for the bouncing ball contribution to the 
fluctuating part of the staircase function for an electro- 
magnetic cavity in three dimensions: 

^et' flUC W = S( E ~ n ") ~ I*" + \ X2 ) 

(A9) 

with X = kl/n. 
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FIG. 1. Geometry of the 3D-cup. The investigated mi- 
crowave resonator together with its measures (in mm) is given 
on the left side. The model of the barrel billiard which was 
used for several numerical simulations is shown on the right 
side. The latter has been constructed from one half of a 
three-axial ellipsoid cut with an additional plane at distance 
h from the origin. 

FIG. 2. Measured transmission spectrum in a range be- 
tween 15 and 20GHz. The ordinate shows the ratio of the 
output power relative to the input power on a logarithmic 
scale. As can be seen from the figure the eigenfrequencies 
appear as sharp peaks with a Q-value of up to 10 5 and a 
signal-to-noise ratio of up to 50dB. 

FIG. 3. Remaining fluctuating part of the stair- 
case function after the extraction of the smooth part 
N amooth (f) = Vi ■ f + V 2 ■ f + V 3 . The particular form of the 
oscillations around zero indicates that the smooth part was 
described consistently. 



FIG. 6. Poincare surface of section for the classical barrel 
billiard. The figure shows the resulting patterns for the conju- 
gated variables (z,p z ) after 16000 collisions with the bound- 
ary. Note that the momentum \p\ is normalized to unity. As 
can be seen the phase space is split in regular stripes and a 
chaotic sea. The former are produced by a special class of 
stable orbits one of which is shown in projections in the lower 
box of the figure producing one certain strip pointed by an 
arrow. These orbits do not hit the bottom of the cup, thus 
they only "see" the regular ellipsoid. 

FIG. 7. Length spectrum (upper part) from the Fourier 
transformed fluctuating part of the eigenmode density of the 
3D-cup (solid curve) and some numerically simulated peri- 
odic orbits for the barrel billiard (lower part). The lengths 
of these orbits and of their multiples are indicated by arrows 
below the abscissa for comparison with the experiment. In 
addition the contribution of the 3D-bouncing ball orbit with 
hbo = 0.1122m has been extracted from the spectrum (dashed 
curve), the larger arrows above the curve indicate multiples 
of l bbo to guide the eye. 



FIG. 4. Ordinary (upper part) and cumulative (lower part) 
nearest neighbour spacing distribution for the set of 286 mea- 
sured eigenfrequencies. Beside the experimental histograms 
also the curves for pure Poissonian as well as pure GOE char- 
acteristics are shown. 

FIG. 5. Number variance E 2 and spectral rigidity A3 for 
the measured spectrum (circles). Again the limiting curves 
for pure Poissonian and Gaussian distributions are given in 
the figure. Up to a certain value L m „ the E 2 -curve is well de- 
scribed by a mixing-parameter q « 0.30±S;|o (dashed curve). 
Above L max « 10 the curve shows a clear deviation which is 
theoretically expected because of the finite set of eigenmodes 
and the finite lengths of the shortest periodic orbits. The ex- 
perimental spectral rigidity A3 is only within its error range 
compatible with the mixing-parameter of E 2 and possesses a 
four times larger value of L max as expected [^BJ. For a clear 
representation only every second error bar is shown in the 
figure. 



FIG. 8. Fluctuating part of the staircase function for a reg- 
ular box using the first 20000 eigenfrequencies. In the upper 
part an impressive oscillation, with period A/ « 750MHz, 
which is due to the contribution of the shortest 3D-bouncing 
ball orbit (shown in the insert) can be observed. The the- 
oretical investigation, Eq.(pl|), resulted in the smooth curve 
also superimposed in the upper part of the figure. After the 
extraction of this bouncing ball contribution the oscillations 
around zero revealed a lower magnitude and the correspond- 
ing period A/ vanished from the fluctuations (lower part). 

FIG. 9. Length spectrum for the regular box in a range 
between and 1.6m (solid curve). In this presentation the 
discussed bouncing ball orbit leads to an impressive peak at 
hbo = 0.4m. The positions of its multiples are denoted by 
arrows in the figure. Again the contribution of this orbit 
has been extracted and now the remaining spectrum (dashed 
curve) can be seen to be free of any remnant. 

FIG. 10. Nearest neighbour spacing distribution for the 
regular box. Because of the systematic degeneration of al- 
most all TE- and TM-modes, the calculation was restricted 
to one polarization (9978 TE-modes). Due to the fact that the 
system is desymmetrized per definition the limiting curve for 
the chaotic case is given by a single GOE distribution. As ex- 
pected for the regular box the system shows pure Poissonian 
characteristics. 
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